Variations of Casimir energy from a superconducting transition 
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We consider a five-layer Casimir cavity, including a thin superconducting film. We show that 
when the cavity is cooled below the critical temperature for the onset of superconductivity, the 
sharp variation (in the microwave region) of the reflection coefficient of the film produces a variation 
in the value of the Casimir energy. Even though the relative variation in the Casimir energy is very 
small, its magnitude can be comparable to the condensation energy of the superconducting film, 
and thus causes a significant increase in the value of the critical magnetic field, required to destroy 
the superconductivity of the film. The proposed scheme might also help clarifying the current 
controversy about the magnitude of the contribution to Casimir free energy from the TE zero mode, 
as we find that alternative treatments of this mode strongly affect the shift of critical field. 
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I. INTRODUCTION 

In recent years, new and exciting advances in experi- 
mental techniques |fj prompted a great revival of interest 
in the Casimir effect, over fifty years after its theoretical 
discovery Q (for a recent review on both theoretical and 
experimental aspects of the Casimir effect, see Refs. Q- 
0). As is well known, this phenomenon is a manifesta- 
tion of the zero-point fluctuations of the electromagnetic 
field: it is a purely quantum effect and it constitutes one 
of the rare instances of quantum phenomena on a macro- 
scopic scale. 

In his famous paper, Casimir evaluated the force be- 
tween two parallel, electrically neutral, perfectly reflect- 
ing plane mirrors, placed a distance L apart, and found 
it to be attractive and of a magnitude equal to: 



F c = 



hcTT 2 A 

240 L 4 



(1) 



Here, A is the area of the mirrors, which is supposed 
to be much larger than 1? , so that edge effects become 
negligible. The associated energy Ec 



E< 



hen 2 A 
720 L 3 



(2) 



can be interpreted as representing the shift in the zero- 
point energy of the electromagnetic field, between the 
mirrors, when they are adiabatically moved towards each 
other starting from an infinite distance. The Casimir 
force is indeed the dominant interaction between neutral 
bodies at the micrometer or submicrometer scales, and 
by modern experimental techniques it has now been mea- 
sured with an accuracy of a few percent (see Refs. p] and 
Refs. therein). 

The issue of the energy of vacuum fluctuations, a man- 
ifestation of which is provided by the Casimir effect, is 
of the outmost importance in other areas of Physics, 
most notably in Cosmology, where it is known to raise a 
number of problems, when its gravitational effects are 



considered Q. For this reason, it is clearly very im- 
portant to devise laboratory tests to study the physical 
consequences of vacuum fluctuations, as a way to check 
our theoretical understanding of this intriguing quantum 
phenomenon. 

All experiments on the Casimir effect performed so far 
measured the Casimir force, in a number of different ge- 
ometric configurations. In a recent Letter Q we have 
found that by realizing a rigid cavity, including a super- 
conducting film, it should be possible to measure directly 
the variation of Casimir energy accompanying the tran- 
sition to the normal state of the superconducting film, in 
an external magnetic field. 

Apart from the fundamental interest of a direct mea- 
surement of a (variation of) Casimir energy, rather than 
a force, our scheme has the further advantage that, be- 
ing based on rigid cavities, it should make it easier to 
study the dependence of the Casimir effect on the geom- 
etry of the cavity. Indeed, since this effect arises from 
long-range correlations between the dipole moments of 
the atoms forming the walls of the cavity, that are in- 
duced by coupling with the fluctuating electromagnetic 
field, the Casimir energy depends in general on the geo- 
metric features of the cavity. For example, we see from 
Eq. J5J that, in the simple case of two parallel slabs, the 
Casimir energy Ec is negative and is not proportional to 
the volume of the cavity, as would be the case for an ex- 
tensive quantity, but actually depends separately on the 
area and distance of the slabs. Indeed, the dependence 
of Ec on the geometry of the cavity can reach the point 
where it turns from negative to positive, leading to re- 
pulsive forces on the walls. For example (see Ref. 0), 
in the case of a cavity with the shape of a parallelepiped, 
the sign of Ec depends on the ratios among the sides, 
while in the case of a sphere it has long been thought to 
be positive. It is difficult to give a simple intuitive expla- 
nation of these shape effects, as they hinge on a delicate 
process of renormalization, in which the finite final value 
of the Casimir energy is typically expressed as a differ- 
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ence among infinite positive quantities. In fact, there 
exists a debate, in the current literature, whether some 
of these results are true or false, being artifacts resulting 
from an oversimplification in the treatment of the walls 
|10|. Under such circumstances, we think it would be de- 
sirable to have the possibility of an experimental check 
of these statements, for the maximum possible variety of 
configurations. 

Another reason of interest in our superconducting cav- 
ities, not discussed in Q, is that their study will help 
clarifying one of the most debated issues in the current 
literature on the Casimir effect, i.e. the question of the 
contribution from the TE zero mode to thermal correc- 
tions in real metals, i.e. characterized by a finite conduc- 
tivity. At the present time, there is no general agreement 
between the experts on how to deal with this problem, 
and we address the reader to Refs. ^Jj for a survey of 
the extended literature on this topic. A superconduct- 
ing cavity is a very good tool to explore this problem, 
because the variation of Casimir energy across the tran- 
sition precisely arises from the fact that the film has a 
finite conductivity, that changes when it becomes super- 
conducting. Indeed, the computations presented in Sec. 
IV show that the critical field required to destroy the 
superconductivity of our cavity is very sensitive to the 
contribution of the TE zero mode, at the level that the al- 
ternative treatments proposed in the literature may lead 
to predictions that can differ even by 100 percent. 

The plan of the paper is as follows. In Sect. II we 
present the general scheme of superconducting Casimir 
cavities, in Sec. Ill we discuss the applicability of Lif- 
shitz theory [T3 | to the computation of the variation of 
Casimir energy in the superconducting transition. In Sec. 
IV we present the results of our numerical computations, 
including in Subsec. C a study of the contribution from 
the TE zero mode and finally, in Sec. V, we draw some 
conclusions and outline possible future work. We have in- 
cluded two Appendices to review the basic facts on type I 
superconductors that are useful to understand the prop- 
erties of our superconducting cavities. 



II. THE SUPERCONDUCTING CAVITY 

The Casimir cavity that we consider is the planar five- 
layer system depicted in Fig. 1: a thin superconducting 
film of thickness D is placed between two thick metallic 
slabs (made of a non-magnetic and non superconduct- 
ing metal), that constitute the plates of the cavity. The 
gaps of width L separating the film from the plates are 
supposed to be filled with some insulating material. 

We consider cooling the cavity at a temperature T be- 
low the critical temperature T c for the onset of supercon- 
ductivity in the thin film. A magnetic field parallel to the 
film is then turned on, and its intensity is gradually in- 
creased, until it reaches the critical value H c » , for which 
superconductivity of the film is destroyed. The question 
we ask is whether the Casimir energy stored in the cavity 
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FIG. 1: Scheme of the superconducting five-layer cavity. 



affects in an appreciable way the value of ff c || . To ad- 
dress this question (see Appendix I for a short review of 
the magnetic properties of type I superconductors), we 
start from Eq. (|A10I) that connects H c « to the (density 
of) condensation energy e(T) of the film. Upon multiply- 
ing both sides of Eq. (|A10|) by the volume V occupied 
by the film, we obtain: 



V_ / g e ||(T) 
8n { p 



— ^cond(^) 



(3) 



where £ CO nd(T) = V e(T) is the condensation energy of 
the film. When the film is included in a cavity, the r.h.s. 
of the above Equation should be augmented by an extra 
term, i.e. the AFg (T) representing the difference 

AFP(T):=FP(T)-FP(T) (4) 

among the Casimir free energies of the cavity, in the su- 
perconducting (s) and in the normal (n) states of the film 
(in zero magnetic field). While we postpone to the next 

Section the computation of AF^ (T), we observe here 
that the presence of such a term should be expected on 
physical ground, because the transition to superconduc- 
tivity causes a sharp variation in the reflective power of 
a film in the microwave region of the spectrum |13| , and 
this obviously produces a change in the Casimir energy, 
which depends on the reflection coefficients of the plates. 
Thus, for the cavity case, Eq. should be replaced by 



£ cond (T) + AF { E C) (T) . (5) 



In writing this equation, the assumption has been made 
that the value of the coefficient p does not change when 
a film is placed in a cavity. This approximation amounts 
to saying that all quantities referring to the film, like the 
penetration depth, condensation energy etc. are not af- 
fected by virtual photons in the surrounding cavity. This 
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is a very good approximation, since the leading effect 
of radiative corrections is a small renormalization of the 
electron mass 0] of order e := a x huj c /(mc 2 ) (up to 
logarithmic corrections), where a is the fine structure 
constant, lu c = c/L is the typical angular frequency of 
virtual photons, and m is the electron mass. Even for L 
as small as 10 nm, e ~ 3 x 10~ 7 . The associated shift of 
critical field is of the same order of magnitude, and thus 
negligible with respect to that caused by AF^\ which 
will turn out to be of some percent (see Sec. IV). 

Equation JSJ deserves further comments. Consider 
first the relative shift of parallel critical field: 



5H, 



H c 



H, 



c|| 



(6) 



which, according to Eqs. 10 and ©, is equal to 



SH C 



AF, 



(C) 



1 - 1 



(7) 



Now, the relative variation AF^ / Fn of Casimir en- 
ergy should be expected a very small number, because 
the Casimir energy is mostly associated with photon en- 
ergies of order Tiuj c , i.e. 10 -j- 20 eV for L ~ 10 nm, while 
the transition to superconductivity affects the reflective 
power only at photon energies of order kT c w 10 -4 eV 
for T c ~ 1 K. However, and this is the key point, Eq. 
(|7|) shows that the shift of critical field is determined by 

(C) 

the variation of Casimir free energy AF^, as compared 
to the condensation energy £ C ond of the film. The lat- 
ter is very small: for a thin film with, say, an area A of 
1 cm 2 (the area A is not really important, because for 
A ^> L 2 both £ con d and AF^ are proportional to A) 
and a thickness of a few nm it is easily of order 10 -8 erg. 
By contrast, we see from Eq. © that a typical Casimir 
energy for a cavity with an area of 1 cm 2 and a width 
L = 10 nm, has a magnitude of 0.43 erg, i.e. over seven 
orders of magnitude larger than the condensation energy 
of the film. This implies that even a relative variation 
of Casimir energy as small as one part in 10 s would still 
correspond to more than 10% of the condensation energy 
of the film, and would induce a shift of critical field of 
over 5 %. 

The huge increase in sensitivity that we predict by 
virtue of a second energy scale, the condensation en- 
ergy, which is orders of magnitude smaller than the scale 
characteristic of the Casimir effect, is the marking differ- 
ence between our scheme and previous attempts at a di- 
rect measurement of the relative variation of the Casimir 
force, caused by a change in the reflective power of the 
mirrors. For example, Ref. [l5l | reports on an experi- 
ment, based on the technology of hydrogen switchable 
mirrors. Despite the large modulation of reflectivity in 
the optical region of the spectrum, which is very relevant 
for typical submicron separations between the mirrors, no 
detectable variation in the Casimir force could be seen, 



possibly due to lack of the necessary sensitivity in force 
measurements . 

Another comment on Eq. is in order. As we shall 

see in Sec. IV, and what is a distinctive feature of the 

(c) 

Casimir effect, the variation of Casimir energy AF^ 
turns out to depend on L, while the condensation energy 
does not. This implies an L dependence of the shift of 
critical field, which, if seen, would be strong evidence for 
the Casimir origin of the effect. In this respect, we notice 
that having chosen for the outer plates a non supercon- 
ducting material, we do not have to worry about possible 
L-dependent Josephson couplings among the film and the 
lateral plates, that could produce a similar effect. 



III. CALCULATION OF THE VARIATION OF 
CASIMIR ENERGY IN THE 
SUPERCONDUCTING TRANSITION 

In this Section, we show how to compute the differ- 
ence among the Casimir free energies of the cavity, in 
the superconducting (s) and in the normal (n) states of 
the film (in zero magnetic field). In this computation 
we shall use the generalization to multilayer systems of 
the theory developed by Lifshitz |l2j to compute the van 
der Waals and Casimir energy for a plane-parallel cav- 
ity formed by two thick dielectric walls separated by an 
empty gap. 

Before we turn to this computation, it is however neces- 
sary to discuss first the applicability of the Lifshitz theory 
to a superconducting cavity. Another concern is to ascer- 
tain whether we can use that theory to reliably evaluate 

( o 

AF^ (T) for separations L and for film thicknesses D as 
small as a few nm, which are the values of interest for us. 
Now, the main assumption of Lifshitz theory is that, in 
the relevant range of frequencies and wave vectors, one 
can describe the propagation of electromagnetic waves 
inside the media forming the cavity, in terms of a com- 
plex permittivity, depending only on the frequency to. It 
is important to stress that the theory includes also non- 
retarded effects, and hence it has as limiting cases both 
van der Waals forces (that become important at small 
distances, like those we consider) and Casimir forces. On 
this ground, it has been used recently to study van der 
Waals interactions among thin metal films (of thickness 
about 10 A), till very small separations (a few A) [lt| . 

When spatial non-locality effects become important, so 
that the complex permittivity becomes a function also of 
the wave vector q, the Lifshitz theory is no longer appli- 
cable 0](see a l so the discussion on p. 126 in Ref. 0]). 
Now it is clear that non-local effects are important, in 
general, in superconductors, by virtue of the very small 
skin depth (of the order of the penetration depth A, typ- 
ically 500 A or so) also in the frequency region charac- 
teristic of the normal skin effect in normal metals (for an 
interesting discussion of non-local effects in the compu- 
tation of dispersion forces in superconductors, see Ref. 
[ill'). Moreover, for small separations L, as well as for 
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small film thicknesses D < 25 — 30 nm and/or at cryo- 
genic temperatures non-local effects become important 
also in the normal state , and indeed it has been advo- 
cated that space dispersion should be taken into account, 
for example, to evaluate the influence of thin metal coat- 
ings, that are used to protect the plates in most of the 
current experiments on the Casimir effect (see [3j and 
last of Refs. p^j . The influence of space dispersion on 
the Casimir effect at room temperature is studied also in 
the recent paper |2jj). These considerations lead to the 
conclusion that, for a reliable evaluation of the individual 

(O (G) 

Casimir free energies F n (T) and Fr ; (T), it is neces- 
sary to consider space dispersion. However, and this is 
a key point, space dispersion is unimportant for the pur- 
pose of computing the difference AF^, (T) between the 
Casimir energies in the two states of the film, which is the 
only quantity of interest for us. The reason is that the 
optical properties of thin films (with a thickness D much 
smaller than the skin depth or correlation length £), in 
the normal and in the superconducting states, are indis- 
tinguishable for photon energies larger than a few times 
kT c , as accurate measurements have shown 13]. This 
implies that, in the computation of AF^ , the only rel- 
evant photon energies are those below a few times k T c 
(corresponding to the far IR), which is where the optical 
properties of the film actually change when it becomes 
superconducting. In this wavelength region, the quoted 
experiment shows that the transmittivity data for thin 
superconducting films can be well interpreted in terms 
of a complex permittivity that depends only on the fre- 



where A 3> L 2 is the area of the cavity, kj_ = (fci, k 2 ) de- 
notes the two-dimensional wave vector in the xy plane, 
while p T i u> kf S p TE ) denote the proper frequen- 
cies of the TM (TE) modes, in the n/s states of the film, 
respectively. 



where A^} (i() is the expression in Eq. (4.7) of (eval- 
uated for £q = e n / s ) and A^j soo (iQ denotes the asymp- 



quency, and is independent of the film thickness. 

Having established the applicability of the Lifshitz the- 
ory to our superconducting cavity, we recall that in the 
original treatment by Lifshitz, dealing with two thick di- 
electric slabs separated by an empty gap, the Casimir 
force was obtained by evaluating, at points in the gap re- 
gion, the Maxwell stress tensor associated with the elec- 
tromagnetic fields generated by the randomly fluctuating 
currents in the interior of the dielectric slabs. Although 
physically transparent, this approach has never been ex- 
tended to multilayer systems, due to the complexity of 
the computations that are involved. The desired gener- 
alization can however be obtained by using an alterna- 
tive approach, similar to the one used by Casimir him- 
self in his study of an ideal metallic cavity, in which the 
Casimir energy is obtained as the zero-point energy of 
the electromagnetic field inside the cavity. This method 
was applied to dielectric cavities in Refs. |2l|, where it 
was shown that one recovers the Lifshitz result , an d was 
later generalized to multilayer systems in Ref. |22l |. 

We thus consider the five-layer system depicted in Fig. 
1. The electric permittivities of the layers are denoted 
as follows: e n i s represents the permittivity of the film, in 
the n/s states respectively, while t\ is the permittivity of 
the insulating layers. Last, €2 is the permittivity of the 
outermost thick normal metallic plates. 

We consider first the T — case. Then, the Casimir 
free energy coincides with the Casimir energy E^ , and 
we can write the unrenormalized variation of Casimir en- 
ergy AE { C) (L,D) as 



(8) 

I 

Upon using the argument theorem, and by subtracting 
the contribution corresponding to infinite separation L 
(for details, we address the reader to Chap. 4 in Ref. [3j) 
[35J , we can rewrite the renormalized sums in Eq. J3J) as 
integrals over complex frequencies i(,\ 



(9) 

I 

totic value of A^j s (iQ in the limit L — > 00 (correspond- 
ing to the limit d — * 00 with the notation of Q)- A sim- 
ilar expression can be written for the TE modes, which 



I 

AP (C),r n >, _ 4 h f dk 1 dk 2 (n,TM) ,(n,TEU -ST^f (s,TM) (s,TB)A 



\ ^ (n, TM) \ ^ (s,TM) 



1 f°° „ (y A n\iC) , Ai 1} «) 
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involves the quantity A^l(iQ denned in Eq. (4.9) of 0- sion for the (renormalized) variation AE^ C \L, D) of the 
Upon inserting Eq. © , and the analogous expression for Casimir energy: 
TE modes, into Eq. JBJ one gets the following expres- 



AE^ = A 



2tt 



-.TE 



log 



Q 



TE 



log 



Q 



TM 



Q 



TM 



(10) 



where we set 



Qi 



(TM/TE) 



(0 



= Ai 1/2) (z0 

(<0 



A 



(1/2) 



I = n,s . 



(11) 



Upon performing the change of variables k\ — (p 2 — 
l)£ 2 /c 2 in the integral over k±, the above expression for 
AE( C \L,D) turns into 



AE^ = 



%A 

47T 2 C 2 Ji 



f)TE n TM 

pdp I dCC 2 (logS- + log^ , 



Q 



~)TM 



(12) 



where the explicit expression of the coefficients q(™/ tb ) j s 



^ _ (1 _ A TE,TM^TE,TM e _ KKlL/c)2 ^TE/TM _ tfE/TM^ K L/c^K, D/c 



1 - {A\f l ™fe- 2C - K ' D I C 



TE K.j K, 
01 ' A 7 +A, 



j A TM _ K 3 e l «) - K l e 3 «) 



K, e, (»C) + K t ej «) 



A', 



^e,^C)-l+P 2 , I = n,s ; j, / = 1,2, n, s. (14) 



Formulae along similar lines, involving the reflection 
coefficients of the layers, have been obtained for example 
in Refs.[13. 

The extension of the above formulae to the case of finite 
temperature is straightforward. As is well known this 
amounts to the replacement in Eq. (|10|) of the integration 



J d(/2ir by the summation kT/HJ^i over ^ ne Matsubara 
frequencies £/ = 2irl/p, where (3 = h/(kT), which leads 
to the following expression for the variation AF E (T) of 
Casimir free energy: 



J 



AF { E C) {T)=A 



kT_ 

2 ^ 

1— — O0 



log 



Q 



TE 



-)TE 



log 



Q 



TM 



Q 



TM 



(15) 



r 



As we see, Eqs. I|12ll5ll involve the electric permittivities 
e (iQ of the various layers at imaginary frequencies z£. 
For these functions, we have made the following choices. 



For the outermost metal plates, we use a Drude model 



for the electric permittivity: 
e D (uj) = 1 



O 2 



where f2 is the plasma frequency and 7 



(16) 
1/r, with 
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r the relaxation time. We denote by and T2 the 
values of these quantities for the outer plates. As is 
well known, the Drude model provides a very good ap- 
proximation in the low-frequency range uj w 2k T c /h ~ 
10 11 -j- 10 12 rad/sec which is involved in the computation 
of AF^ C) (T) and AE^ . The relaxation time is tem- 
perature dependent and for an ideal metal it becomes 
infinite at T = 0. However, in real metals, the relax- 
ation time stops increasing at sufficiently low tempera- 
tures (typically of the order of a few K), where it reaches 
a saturation value, which is determined by the impurities 
that are present in the metal. Since in a superconducting 
cavity the temperatures are very low, we can assume that 
T2 has reached its saturation value and therefore we can 
treat it as a constant. The continuation of Eq. l|ltil) to 
the imaginary axis is of course straightforward and gives 



e D (iQ = 1 + 



n 2 



C(C + 7) ' 



(17) 



For the insulating layers, we take a constant permit- 
tivity, equal to the static value: 



ei(w) = ei(0) . 



(18) 



Again, this is a good approximation in the range of fre- 
quencies that we consider. 

As for the film, in the normal state we use again the 
Drude expression Eq. I|16|l . with appropriate values for 
the plasma frequency Q n and the relaxation time r n . 

The permittivity e s (iC) of the film in the superconduct- 
ing state cannot be given in closed form and we have eval- 
uated it as explained in Appendix B, starting from the 
Mattis-Bardeen formula for the conductivity ct s {oj) of a 
superconductor in the local limit q — > of BCS theory. 



IV. NUMERICAL COMPUTATIONS 

In this Section we present the results of our numeri- 
cal computations. For the convenience of the reader, the 
material is divided into four separate Subsections. In the 
first Subsection we analyze the qualitative features to be 
expected from a numerical evaluation of Eqs. i|10|) - (|15[) . 

In Subsec. B we present the computation of the varia- 

(c) 

tion AFjf' of Cas imir free energy in the superconducting 
transition. In Subsect. C we discuss the controversial is- 
sue of the contribution arising from the TE zero mode. 
Finally, in the last Subsection, we present our results for 
the shift of critical magnetic field. 



A. Qualitative features of AF^ 

In this Subsection we discuss some qualitative features 
of Eqs. IjlUp 1151) that are useful as a check of the nu- 
merical computations. We recall first of all that, with our 
choices for the cavity geometry and for the permittivities 



(c) 

of the various layers, the value of AF^ depends on the 
following sets of parameters: 

1. the thickness of the film D and the width of the 
insulating gaps L (see Fig. 1); 

2. the temperature T, which we shall express in units 
of T c , by means of t = T/T c ; 

3. the plasma frequency Q2 and the saturation value 
T2 for the relaxation time of the outer plates; 

4. the dielectric constant ei(0) of the insulating gaps; 

5. the plasma frequency f2 n and the saturation value 
T n for the relaxation time of the film, in its normal 
state; 

6. the BCS gap A(T) of the superconducting film. 

Some simple properties of our expression for AF^ (T) 
are obvious. First, it is clear from Eq. (|13|l that 



or /TM (c, P ) 



1 for L — > 00, and thus AF 







in this limit, as it should. Another important limit is 
that for L — ► 0. In this limit, AF^ approaches a finite 
value, because, as pointed earlier, for large C the permit- 
tivities of the film in the normal and in the superconduct- 
ing states become undistinguishable, and then the ratios 
qTE/tm jqTE/tm a pp roacn one ( anc [ their logarithms 

approach zero) faster than Q^f s , and this ensures 
convergence of Eqs. i|12|) and i|15[l also for L = 0. 

We consider now the dependence of Ai^ C) (T) on D. 
For D — > 00, the expressions of Q^ E ^™ reduce to 

lim Q^ 1 ™ = (1 - A ^/ TM Af 2 £/TM e- 2 ^ L ^f , 

D—*oo 

which concides with the square of the analogous expres- 
sion for a three-layer cavity consisting only of the outer 
plate, the insulating layer and an infinite superconduct- 
ing layer. Upon taking the logarithm of Q^ E ^™ we 
(c) 

thus obtain that AF^ is twice the value for the three- 
layer system, which is what one would have expected 
since in the limit D — > 00 our five-layered cavity be- 
comes equivalent to a system of two identical three- 
layered cavities. The rate of approach to the limiting 
double three-layer system is controlled by the exponential 
factors exp(— 2DK n / s ( ! /c) 1 which set the scale Do above 
(c) 

which AFg becomes independent of D and reaches its 
D — > 00 limiting value. The value of Dq can be estimated 

as follows. For photon frequencies C of order £ = 2A/h, 

(c) 

which are the relevant ones for AF^, , it turns out that 
e s (i£) > e n (iQ (see Appendix B), and therefore we have 
that K s (C,p) > K n (C,p) > K n ((, 0). On the other hand, 
using the Drude model Eq. ljT7Jl for e n {K) we estimate 
(K n ((, 0) w C \/e««) ~ + y, which implies 



cxp(-2DK n/s (/c) < exp 



2Dn n x 
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We see from this Equation that the exponential becomes 
negligible for D ^S> D , where 



D 



2 ft* 



(20) 



For a typical value of the plasma frequency h£l n ~ (10 -j- 
20) eV, we thus obtain 

Dq w 10 nm . 

The above number is very important for us, because 
it provides an estimate of the thickness D of the film, 
for which the largest shift of critical field should be ex- 
pected. Indeed, we see from Eq. (0 that the largest 
shift is obtained for the value of D that maximizes 
the ratio AFg /£ CO nd- Now, the condensation energy 
£cond = V e(T) is proportional to the volume V of the 
film and then to D (in reality, for certain superconduc- 
tors and for ultrathin films with D « the critical 
temperature T c and then the density of condensation en- 
ergy e(t) may depend on D. See discussion on Be films 
in the next Subsection). Therefore, for D ^> Dq, when 

AF^ becomes independent of D, the ratio AF^ /£ con< \ 
decreases like 1/D because of the unlimited increase of 
£cond with D. It then follows that the optimal values of 
D must be of order Dq or so, i.e. of a few nm, and this 
is the range that we shall consider in what follows. 



B. Numerical computation of the variation of 
Casimir free energy 

The variation of Casimir free energy AFg(T), Eq. 
(JT5J, has been computed numerically. For a supercon- 
ducting film with a critical temperature T c , the gap A(T) 
has been estimated by using the approximate BCS for- 
mula Eq. (|A13(1 . In what follows we present some of 
our results in terms of the so-called impurity parameter 
y — h/ (2A(T)t„), which is commonly used, in place of 
the relaxation time r„, as a convenient measure of the 
degree of purity of the superconductor. In the range of 
parameters that we consider, the value of AF^ is in- 
dependent, to better than four significant digits, of the 
value of ei(0), for 1 < ei(0) < 300. In all computations 
presented below, we fix once and for all ei(0) = 100. 

Another important remark is that, in all cases that we 
considered, the contribution of TM modes to AFg is 
completely negligible, being three or four orders of mag- 
nitude smaller than that of TE modes. We postpone to 
next Section a discussion of the implications of this fact. 

In Fig. 2 we show the plot of AF^ (in erg) as a 
function of the width L (in nm) of the insulating gap, for 
D = 5 nm, T c = 0.5 K, y = 15, t = 0.9, ft n = ft 2 = 18.9 
eV, r 2 = 2.4 x 10~ 12 sec. We observe that AF^ ] is al- 
ways positive, which corresponds to the intuitive expecta- 
tion that transition to superconductivity of the film leads 
to a stronger Casimir effect, i.e. to lower Casimir free en- 
ergy. The data can be fit very accurately by a curve of 
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FIG. 2: Plots of AF^I 



in erg) as a function of L (in nm) 
for D — 5 nm, T c — 0.5 K, t = 0.9. See text for the values 
of the other parameters. Also shown is the plot (dashed line) 
of a fit of the type 1/(1 + (L/L ) a ), with L = 8.3 nm and 
a = 1.15. 



the type 



a4 c) (io« 



i 



1 + (L/L ) a 



(21) 



where Lq = 8.3 nm and a — 1.15. We thus see, as an- 
ticipated earlier, that AF^ (L) approaches a finite limit 
for L — > 0. This feature of AF^ is helpful, from an 
experimental point of view, because it implies that the 

possible roughness of the film and of the lateral plates will 

(C) 

have little influence on the value of AF^ , for L suffi- 
ciently small. By contrast, we recall that the roughness 
of the surfaces is one of the main sources of uncertainty 
in the theoretical analysis of standard measurements of 
the Casimir force, for submicron separations Q . 

In Fig. 3 we show the plots of AF^' (in erg), as a 
function of the impurity parameter y, for the following 
three combinations of values of the plasma frequencies 
for the film and the lateral plates: ft„ = ft 2 = 18.9 
eV (solid curve), ft„ = 18.9, ft 2 = 12 eV (point-dashed 
curve) and ft„ = ft 2 = 12 eV (dotted curve) (the val- 
ues of 18.9 eV and 12 eV correspond respectively, to Be 
and Al). All curves have been computed for L = 10 nm, 
D = 5 nm, T c = 0.5 K, t = 0.9, r 2 = 2.4 x 10" 12 sec. We 
see that AF^ has a maximum for y ~ 10 -j- 15. While 

the decrease of AF^ for large values of y is obvious, it 
is perhaps surprising to see that it decreases also when 
y becomes very small, i.e. when the degree of purity of 
the film is improved. In fact, this latter behavior is ex- 
plained by observing that, for small values of y, the film 
is so pure that already in the normal state, it behaves as a 
very good conductor and then transition to superconduc- 
tivity cannot increase significantly the Casimir energy. 
(c) 

Notice also that AF% is very sensitive to the values for 
the plasma frequencies both of the film and of the lat- 
eral plates, and it increases sensibly when either one is 
increased. 

(c) 

In Fig. 4 we show the plot of AF^ as a function of 
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FIG. 3: Plots of AF { E C) (in erg) for L = 10 nm, D = 5 
nm, T c — 0.5 K and t = 0.9, as a function of the impurity 
parameter y, for O n = fl z = 18.9 eV (solid line), fl n = 18.9, 
Q2 = 12 eV (dashed line) and for fi n — fl^ = 12 eV (point- 
dashed line). See text for values of the other parameters. 



FIG. 5: Plot (solid line) of AF E C ' (in erg) as a function of 
t for L = 10 nm, D = 5 nm, T c = 0.5 K, r n = 5 x 10" 13 
sec. The point-dashed line includes the contribution of the 
TE zero mode, computed using the plasma model for the 
lateral plates. Also shown (dashed line) is the plot of the low- 
temperature limit of the Matsubara sum, Eq. 1121 . See text 
for further details. 
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FIG. 4: Plots of AF ( z } (in erg) as a function of T c (in K) 
for L = 10 nm, D = 5 nm, y = 15 and t = 0.9. See text for 
values of the other parameters. 



the critical temperature T c of the film, for L = 10 nm, 
D=5 nm, and t = 0.9. All other parameters are as in 
Fig. 



2. The figure shows that AFg depends linearly on 



In Fig. 5 we show the plot of A-Fg as a function of 
t, for T c = 0.5 and r„ = 5 x 10~ 13 sec (all other pa- 
rameters are as in Fig. 4). Notice that in this plot we 
are holding r„ constant and hence, by virtue of the tem- 
perature dependence of the gap A(T) (see Eq. (|A13ll ). 
the impurity parameter y is not constant (for T = 0, we 

have y(0) = 8.7). Note also that, for T/T c -> 1, AF { f } 
approaches zero linearly in (1 — t). Also shown in the 
Figure (dashed line) is the plot of the low-temperature 
approximation to AF^ , obtained by replacing the Mat- 
subara sum Eq. I|15fl by the integral over (, as in Eq. 
(|12fl . This approximation was used in Ref. 0, and as 
we see from Fig. 5 it works fairly well in the whole range 
of temperatures. For the explanation of the third curve 
(point-dashed line) shown in Fig. 5, see next Subsection. 



C. Contribution from the TE zero mode. 

In the current literature on thermal corrections to the 
Casimir effect there is an ongoing controversy concerning 
the proper way of calculating the contribution of the TE 
zero mode (i.e. the I — term in the Matsubara sum) 
to the Casimir free energy (see Refs. for a detailed 
discussion of different points of view on this problem). 
The essence of the question is whether one can use the 
Lifshitz theory to evaluate the TE zero mode, and in the 
affirmative case, for what choice of the metal permittivity 
function e(i£). Indeed, all the trouble arises from the fact 
that the computation of this mode involves the quantity 
C: 



C:=lim(C 2 e(iC)) 



(22) 



If one uses the Drude function Eq. I|17fl (with a finite 
value of r) one obtains C = 0, and this implies that the 
TE zero mode gives zero contribution to the Casimir free 
energy, irrespective of how large r is. On the contrary, if 
one uses the simpler plasma model that gives 



(iC) 



Q 2 

p 

c 2 



(23) 



one finds C = fi 2 ,, and then the zero mode gives a non 
vanishing contribution, reproducing the ideal metal case 
in the limit of infinite plasma frequency. While we ad- 
dress the reader to Refs. ^lj for a discussion of the rea- 
sons in favor of the alternative approaches to this ques- 
tion that have been proposed by a number of authors, we 

would like to examine the consequences of this problem 

( c) 

for our computation of AF^ . 

All computations presented in the previous Section 
used the Drude model for the permittivities for both the 
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lateral plates and the film, in the normal state. Since, 
as pointed out earlier, the contribution from the TM 

(c) 

modes is negligible, it follows that the values of AFg 
that we computed receive no contribution from the TE 
zero mode. 

It is interesting to see what happens in our computa- 
tions if, in the evaluation of the zero mode, we replace 
the Drude model by the plasma model. We shall do it 
only for the lateral plates. For the superconducting film 
we prefer keep using the Drude model, since it is the one 
to which the BCS formula of the permittivity converges 
for T -> T c . Thus, if we denote by AF [ p{TE z.m.) the 
contribution of the TE zero mode, we obtain for the total 
variation of Casimir free energy the expression 



AF 



(C) 



AF^\TE z.m.) + AF 



(C) 



(24) 
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FIG. 6: Plot 
function of t. 



of the ratio AF^ C) (TE z.m.)/AF^ c) 



outer plates. 



All parameters are same as in 



as a 
Fig. 5. 



AFg'(TE z.m.) is evaluated using the plasma model for the 



— (C) 

In Fig. 5 we plot (dashed line) AF E as a 
function of t, while in Fig. 6 we show the ratio 
AFg (TE z.m.)/AFg , again as a function oft (all pa- 
rameters are as in Fig. 5) . We observe that the weight of 
the zero mode increases as we approach T c . The reason of 
this is easy to understand: the zero mode becomes more 

and more important because a decreasing number of Mat- 

(c) 

subara modes contribute to AF S , as one approaches T c . 
This is so because the quantities (Q^ E /Q^ E )(i() that oc- 
cur in Eq. I|15|l are substantially different from one only 
for complex frequencies £ of order of a few tens of times 
kT c /h. Since the l-th Matsubara mode has a frequency 

equal to 2-7T I k T/h, it is clear that the number of terms 

(c) 

effectively contributing to AFg should be roughly pro- 
portional to T c /T, and hence it is large for T <C T c , but 
becomes small for T comparable to T c . Notice that close 
to T c the contribution of the zero mode is practically of 
the same magnitude as that of AF^ and therefore its 
inclusion doubles the variation of Casimir free energy. 



D. Shift of critical field 



Having estimated the variation of Casimir energy 

(c) 

AFg across the transition, we are now ready to com- 
pute the shift of parallel critical field, using Eq. JSJ. The 
result depends of course on the superconductor, and in 
order to get a feeling of what to expect depending on this 
choice, it is convenient to consider the approximate ex- 
pression for the relative shift Eq. J7J. According to this 
formula, the shift is larger when £ con( j is smaller and/or 
AF^ is larger. Now, we see from Fig. 4 that AF^, 
is almost proportional to T c , while from Eqs. I|A4(1 and 
(|X5l) in Appendix A we see £ coa d{T) oc H%(0). There- 
fore, if we forget for a moment other influential factors 
(like the values for the impurity parameter y, the plasma 
frequency and the relaxation time that we shall consider 



TABLE I: Values of T c (in K), H c (0) (in Oe) and s = 
# C 2 (0)/T C in (Oe 2 /K) for some type I superconductors. 





T c 


flc(0) 


H 2 M/T C 


Be 


0.024 


1.08 


48.5 


Zu 


0.88 


53 


3192 


Al 


1.2 


105 


9343 


Sn 


3.72 


306 


25171 



later) we obtain from Eq. Q the estimate 



AF 



(C) 



T, 



5H r 



H c \\ 2£ cond (T) #2(0) 



oc 



(25) 



We thus see that the relative shift should be inversely 
proportional to the quantity s := H^(0)/T c , and so su- 
perconductors with small values of s should be preferable. 
Now, if we use the empirical law Eq. (|A6|) to estimate the 
dependence of H c (0) on the critical temperature of the 
superconductor, we find that s goes like T c 16 , and this im- 
plies that we should first consider superconductors with 
low T c . In reality, things are not so simple. Consider 
for example the case of Be. It has a very low T c = 24 
mK and a correspondingly low H c (0) = 1.079 ± 0.002 
Oe [23, giving s = 48.5 Oe 2 /K, which is much smaller 
than the values of s for other superconductors. As a 
comparison, we quote in Table I the values of T c , H c (0) 
and s for a number of type I superconductors. From these 
data it would seem that Be should give a great advantage 
over other superconductors. However, unfortunately, the 
above low value of T c refers to bulk samples, while very 
thin Be films (with thicknesses of a few nm) have a much 
higher T c . For example, the authors of Ref. j2?J report a 
T c of 0.5 K for a thickness of 5 nm. The value of the con- 
densation energy for such a film is not reported in Ref. 
|27j . but one can estimate that £ con d scales with respect 
to the bulk value in the same proportion as the square 
of the critical temperature (see Appendix A) . If we insist 
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FIG. 7: Plot (solid line) of the relative shift of parallel critical 
field for a Be cavity, as a function of T/T c . The point-dashed 
line includes the contribution of the TE zero modes, evaluated 
as in Fig. 5. All parameters are as in Fig. 5. 



in expressing the condensation energy in the form of Eq. 
(|A5|) . this amounts to saying that the thermodynamical 
field of a thin film scales with its critical temperature. 
This gives us an estimate of H c (0) =22.5 Oe for the films 
of Ref. 0, and thus we find s = 1010 Oe 2 /K, which is 
much worse than the bulk value, but still good compared 
with other type I superconductors. 

We consider now the other factors affecting the rela- 
tive shift of critical field, that we have neglected in the 
above discussion. The first important factor is the value 
of the plasma frequency f2. As we see from Fig. 3, larger 

values of £1 both for the film and the lateral plates lead 

(c) 

to significantly larger values for AFg . Now we recall 
that Be has TiVl = 18.9 eV, which is much larger than the 
value for other metals (for example Al, which has Kl = 
12 eV). Therefore it would be advantageous to use this 
material both for the film and the lateral plates. In this 
respect, the circumstance that the critical temperature 
of a thin Be film is much higher than the bulk value can 
be exploited to our advantage, as it allows to use Be both 
for the film and the lateral thick plates, leaving us with 
a wide range of temperatures for which the film is super- 
conducting, while the lateral mirrors are normal, as we 
require. According to Fig. 3, we can gain in this way 
a factor larger than two, with respect to say, a cavity 
made solely of Al. Another reason in favor of Be is re- 
lated to the fact that the critical field H c \\ of thin films is 
not defined in general as sharply as in bulk samples, and 
one rather has a fuzzy critical region of finite width. It is 
clear therefore that preference should be given to materi- 
als for which this indeterminacy is smallest. It turns out 
that Be can be deposited in ultrathin films, characterized 
by an extraordinary degree of uniformity. A consequence 
of this uniformity is that even for thicknesses D of a few 
nm, Be films possess very sharp parallel critical magnetic 
fields ff c ii , that can be measured with a precision of a few 
parts in a thousand • 

Having collected a number of arguments in favor of 
Be, we can now show the results of our computations 



for the shift of critical field. In Fig. 8 we show (solid 
line) the relative shift of the critical parallel field of a Be 
film, as a function of t. Note that, according to Eq. (J7J 
this quantity is independent of the coefficient p, whose 
precise value is in general not known accurately. The 
Figure has been drawn using the same parameters as in 
Fig. 5. We show also (point-dashed curve) the relative 
shift that is obtained after including the contribution of 
TE zero modes, computed in the way explained in the 
previous Subsection. We see that this inclusion produces 
almost a doubling of the shift. 

Note that the shift is positive, meaning that the critical 
field for the film placed in the cavity is larger than the 
critical field for an identical film outside the cavity. The 

increase of critical field relative shift as one approaches 

(c) 

the critical temperature arises because, for t — > 1, AFg 
and £ con d approach zero at different rates. Indeed, while 
£ C ond vanishes as (1 — t) 2 (see Eq. (|A4H ). we see from Fig. 
5 that AF% vanishes essentially like the first power of 
1 - t. 



V. CONCLUSIONS AND DISCUSSION 

We have studied a superconducting Casimir cavity, in- 
cluding a thin superconducting film and we have shown 
that the variation of Casimir free energy, accompanying 
the superconducting transition, determines a measurable 
shift of the critical magnetic field required to destroy the 
film superconductivity. This constitutes a novel approach 
to the study of the Casimir effect, in which the central 
role is played by the Casimir energy itself, rather than 
the Casimir force, as in all experiments performed so far. 

Our scheme, apart from its novelty, presents a number 
of appealing features. The first advantage stems from the 
fact that it uses rigid cavities. As is well known, the ex- 
perimental difficulty of controlling the parallelism among 
macroscopic plane plates with submicron separations led 
the experimenters to consider simpler geometries that do 
not suffer from this problem, like the sphere-plane one, 
which has been adopted in all precision experiments on 
the Casimir force (with the only exception of the experi- 
ment by Bressi et al. pj , where the plane-parallel config- 
uration is used). This limitation has made it impossible 
so far to explore experimentally one of the most intrigu- 
ing features of the Casimir effect, i.e. its dependence on 
the geometry of the cavity. The use of rigid cavities, as 
in our approach, might make it easier to realize a num- 
ber of geometries that are difficult to realize within the 
standard force-measurements scheme. 

Another interesting feature of our scheme, not men- 
tioned in Ref. [8j, is its sensitivity to the contribution 
from the TE zero mode. This is a controversial issue 
in the current literature on thermal corrections to the 
Casimir effect. It is well known that, for submicron plate 
separations, thermal corrections are negligible at cryo- 
genic temperature, and become relevant only at room 
temperature. However, this is not the case in our su- 
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perconducting cavities, because the shift of critical field 
is completely determined by the Matsubara modes with 
frequencies below or of order kT c /Ti, which is where the 
reflective properties of a film change when it becomes su- 
perconducting. As a consequence, as we have shown in 
Sec. IV C, different treatments of the TE zero mode lead 
to strongly different predictions for the shift of critical 
field, at the level of doubling the shift close to T c , and 
this opens the way to a possible experimental clarification 
of this delicate problem. 

Our computations show that larger shifts of the crit- 
ical field can be obtained by employing for the lateral 
plates and the superconducting film materials with high 
plasma frequencies, which ensure a high reflection coeffi- 
cient. We found also that it is advantageous to use films 
with a value of the impurity parameter y about 10 4- 15. 
Such values of y require films with rather large free mean 
paths, of order 100 nm or so, which can be difficult to 
reach in ultrathin films with thicknesses of 5 4 10 nm, as 
we require. Finally, since larger relative shifts are found 
near the critical temperature, an effort should be made 
to produce ultrathin films with high degree of purity, in 
order to obtain narrow transition regions, with sharply 
defined critical temperatures. 

It is also the case to note that the variation of Casimir 

(c) 

free energy AF^ approaches a finite limit, for vanishing 
separation L of the film from the outer plates. For sep- 
arations of order 10 nm, this behavior of AFg ensures 
a good degree of robustness of the resulting shift of crit- 
ical field with respect to the roughness of the metal sur- 
faces, which is known to give an important contribution 
in standard Casimir force measurements, for submicron 
separations [3j . Equally helpful in this respect is the fact 
that AF^ is independent, to a high accuracy level, of 
the dielectric function of the insulating layers separating 
the film from the outer plates. 

It would be very interesting to obtain an experimen- 
tal verification of the phenomena described in this pa- 
per. Indeed this is the aim of the ALADIN experiment, 
sponsored by INFN, which is currently under way at the 
Dipartimento di Scienze Fisiche dell'Universita di Napoli 
Federico II. Some details about the experimental setup 
are now helpful. In order to detect the shift in the crit- 
ical field of a thin superconducting film resulting from 
the Casimir effect, we foresee to simultaneously deposit 
on a single substrate a certain number of superconduct- 
ing films of same thickness, covered with a thin insulating 
layer, again of same thickness for all films. Half of these 
films will then be covered with a third thick layer of a 
non-superconducting metal, thus forming a three-layer 
Casimir cavity. Such three-layer systems are easier to 
manufacture than the five-layer cavities discussed in this 
paper, and thus, even if they lead to a smaller effect, will 
be considered first. Passage to the five-layer case is ex- 
pected at a later stage, and will be dealt with in a similar 
way. The simultaneous deposition as above is aimed at 
ensuring that all superconducting films, both the single 
ones and those in the cavities, have the same properties 



(thickness, degree of purity etc.). We plan to repeat the 
measurements for various film and oxide thicknesses, in 
order to test our theoretical predictions. 

Besides the authors of the present paper, the exper- 
iment ALADIN has involved D. Born, A. Cassinese, F. 
Chiarella, L. Milano, O. Scaldaferri, F. Tafuri, and R. 
Vaglio. 
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VI. APPENDICES 

For the convenience of the reader, we briefly review in 
these Appendices the basic properties of type I super- 
conductors, that are important for the superconducting 
Casimir cavities to be considered later. Fore a more com- 
plete discussion, we address the reader to the monograph 



APPENDIX A: MAGNETIC PROPERTIES OF 
TYPE I SUPERCONDUCTORS 

As is well known, superconductors tend to expel mag- 
netic fields from their interior (Meissner effect). The re- 
versibility of this effect implies the existence of a well 
defined critical magnetic field H c that destroys supercon- 
ductivity, via a first-order phase transition. If we denote 
by e(T) the difference between the zero- field Hclmoltz 
free energies (per unit volume) f n (T) and f s (T) in the 
normal (n) and the superconducting (s) states respec- 
tively, i.e. 



e(T) := f n (T) - f s (T) , 



(Al) 



(e(T) is called the condensation energy of the supercon- 
ductor) it can be shown that the value of H c is related 
to e(T) by the relation 



H 2 C {T) 
8tt 



e(T) 



(A2) 



Upon using the well known empirical parabolic law for 



H C {T) » H c (0) (1 - t 2 



(A3) 



where t = T/T c with T c the critical temperature, it fol- 
lows from Eq. i|A2|l that the condensation energy has the 
following dependence on t: 



e(T) = e(Q)(l-ty, 



(A4) 
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where 



e(0) 



g c 2 (Q) 

8tt 



(A5) 



Another empirical law of interest |25j is the one giving 
the approximate dependence of H c (0) on T c (valid for 
type I superconductors) 

H c (0) oc T c 13 , (A6) 



and then by using Eq. (|A5|) we obtain an estimate of 
the dependence of e(0) on the critical temperature of the 
superconductor in the form 



e(0) cx T ( 



2.6 



(A7) 



In the case of a thin film, of thickness DcA,( (with 
A the penetration depth and £ the correlation length), 
placed in a parallel magnetic field, expulsion of the mag- 
netic field is incomplete, and consequently the critical 
field increases from H c to H c ii . Unlike the bulk case, the 
transition to the normal phase is of second order, since 
sufficiently thin films and use of the Ginzburg-Landau 
theory show that H c u is proportional to H r (T) 



H 4 (T)=p(T)H c (T). 
The coefficient p has the approximate expression 

9D 2 



P 



7T 6 C 2 



(A8) 



(A9) 



where the second term inside the brackets accounts for 
surface nucleation. Upon using Eq. l|A8p to express H c 



in terms of H c u , we can recast Eq. (IA2|) in the form 



on the preparation procedure. For example, in the case 
of thin Beryllium films one can have a T c as high as 9 
K, to be contrasted with the bulk value of 24 mK [30| . 
In such cases, the condensation energy of the film should 
not be expected to be the same as for the bulk material. 
However, one can estimate it by using the BCS formula 
for the condensation energy (per unit volume) at absolute 
zero: 



e (0) = iiV(0)A 2 (0), 



(A12) 



where N(0) is the density of electronic states of one spin 
in the normal metal at the Fermi surface [3lJ . If we as- 
sume that the gap A(T) is always given by the approxi- 
mate BCS formula yjl 



2A(T) = 3.528 fcT c yi- ^ ^0.9963 + 0.7735^ 

C (A13) 

and that N(0) for the film is not much different from the 
bulk value, we obtain from Eq. I|A12|I 



film 



I film 



bulk 



T} 



bulk 



(A14) 



If we further assume that the temperature dependence of 
the condensation energy of a film is still of the form Eq. 
(|A4|) , the same ratio is found also for the free energies at 
finite values of t. 



APPENDIX B: HIGH-FREQUENCY 
ELECTRODYNAMICS 



1 

8^ 



H C \\{T) 



e(T) 



(A10) 



Recalling that the temperature dependence of A can be 
approximately described by the law 



(All) 



it follows from Eq. I|A10|) that, close to T c , H c \\ ap- 
proaches zero like \/l — T/T c . 

In the case of ultrathin films (with thicknesses of a few 
nm) the magnetic behavior can be quite different because 
the orbital currents are suppressed, and the transition to 
the normal state is driven by coupling of the magnetic 
field to the electronic spins. Under such circumstances, 
the critical field transition can be strongly hysteretic, be- 
low a tricritical point |26j. However, close to T c , there 
is no hysteresis, and if the films are very pure it is still 
possible to observe very sharp transitions, with well de- 
fined critical fields (defined to within a few parts in a 
thousand) |27j . 

It is important to point out that the critical tempera- 
ture T c of a thin film does not coincide in general with 
the bulk value, and it depends on the film thickness and 



As is well known, superconductors show finite dissipa- 
tion when traversed by alternating currents and/or time- 
varying electromagnetic fields. At frequencies uj much 
smaller than the energy gap, Tiu <C A(T), a qualitative 
description of the superconductor response is provided by 
the simple Casimir-Gorter two-fluid model, which leads 
to the following expression for the real part o~'(uj) of the 
complex conductivity o-(lu): 

a' ' (uS) — {i:n s e 2 /2m) S(ui) + (n n e 2 T n /m) (1 +u 2 t 2 )~ 1 , 

(Bl) 

where S(uj) is the Dirac delta function, m and e denote 
the electron mass and charge, n s and n n are the tem- 
perature dependent densities of the superconducting and 
normal electron respectively, and r„ is the residual re- 
laxation time for the normal electrons, as determined by 
the impurities present in the sample. As we see, the nor- 
mal fluid contribution implies a nonzero dissipation at 
all nonzero frequencies. The delta function contribution, 
proportional to n s is a dc contribution, ensuring fulfill- 
ment of the oscillator strength sum rule: 



dto a'(uj) = 



Trne 
2m 



(B2) 



13 



where n = n s + n n is the total electron density. 

An accurate description of the complex conductivity 
of a BCS superconductor, for photon energies up to and 
larger than (twice) the gap, requires that one uses the 
full Mattis-Bardeen theory of electric conduction, which 
is valid also in the domain of the anomalous skin effect 
[3^ |. Taking full account of space dispersion, this theory 
gives rise to an expression of the complex conductivity 
<j s (lo, q) that depends also on the spatial wave vector q. 
In the case of a BCS conductor at a temperature T <T C , 
and in the local limit q — * 0, the resulting expression of 
ct' s (oj) can be written in a form analogous to Eq. I|B1|) : 



a' s (uo) = n5(uj)+a' s (u) 



(B3) 



For lo > 0, (t' s (lo) reads as 



Tine 



poo /> — A 

/ dEJ T + 9{huo - 2A) / dE J D 

J A Ja-Hlu 

(B4)' 



2mujT n 
where 

J T ■= g(uj,T n ,E) 



E + hio E 
tanh — — — tanh 



2kT 



2kT 



(B5) 



Jd ■■= -g(u,T n ,E) tanh I 



with k the Boltzmann constant. Defining 



(B6) 



Pi := y/{E + Hlo) 2 - A 2 , P 2 := V ' E 2 - A 2 , (B7) 



the function g(uj, r„, E) is 



E(E + hiu) + A 2 
E(E + Hlu) + A 2 



1 



J (Pi - P2? + (h/T n ) 2 
1 



PlP 2 J (Pi - P2? + {h/Tnf 



!l = 1 
1 



The coefficient n of the delta function in Eq. 1|B3(I is 
again determined so as to satisfy the sum rule Eq. (|B2|) 
and can be computed exactly according to |3^| 



7TT„A A 

— - — tann 



2kT 



tanh(E/2fcT) 



4A 2 / dE 

I a VE 2 -A 2 [4(E 2 - A 2 ^ 



r 



(Tl/Tn) 2 



(B8) 



The conductivity ct' s (lo) in Eq. I|B3|I can be thought of 
as the sum of three contributions: a 5 function at the 
origin, a broad thermal component that diverges loga- 
rithmically at lo = and a direct absorption compo- 
nent, with an onset at 2A(T) (see the second in Ref. 
[33|]). At any T < T c , complete specification of cr' s (uo) 
requires three parameters: besides the free electron den- 
sity n (or equivalcntly the square of the plasma frequency 
ft 2 = 47rne 2 /m) that provides the overall scale of a' s , and 
the relaxation time for the normal electrons r„, both of 
which already occur in the simple Drude formula Eq. 
(IB9II below, u' s (lo) only depends on one extra parameter, 
i.e. the gap A. In fact, after the frequencies are expressed 
in reduced units x = hu)/(2A), the expression inside the 
brackets on the r.h.s. of Eq. i|B4(l becomes a function 
solely of x, 2A/(fcT) and the so-called impurity param- 
eter y = h/(2AT n ). We point out that this expression 
for a' s is valid for arbitrary relaxation times r„, i.e. for 
arbitrary mean free paths, and in particular it holds in 
the so-called impure limit y = 2A/(hr n ) ^> 1, where the 
effects of non-locality become negligible. We recall that 
the gap A is a temperature dependent quantity. In our 
numerical computations, we used for it the approximate 
formula Eq. I)A13|) . 

We point out that at fixed lo for T — > T c , as well as 
at fixed T < T c for x — > 00, <y' s {u>) approaches the Drude 



1 - • 
0.8 
0.6 
0.4- 
0.2 



1 



FIG. 8: Plots of mo' s {Lu)/(ne 2 T n ), for T/T c = 0.3 (solid 
line), T/T c = 0.9 (dashed line) and T = T c (point-dashed 
line). On the abscissa, the frequency lo is in reduced units 
x = ftw/(2A(0)), and y = 2A(0)/r„ ~ 8.7. 



expression <t' d (uj) 



a' D (co) 



1 



n 2 



Aw 1 + UJ 2 T 2 



(B9) 



The convergence of <j' s (lo) to a' D (uo) in the frequency do- 
main is in fact very fast, and already for x of order 
10 or so a' s becomes undistinguishable from a' D , in ac- 
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cordance with experimental findings [l3| . In Fig. 5, 
we show the plots of cr' s {w) m/(ne 2 T n ), for T/T c = 0.3, 
T/T c — 0.9 and T = T c . The curves are computed for 
yo = 2A(0)/r„ ~ 8.7. Frequencies are measured in re- 
duced units xq = Uw/(2A(0j). 

As we see from Eqs. Ijl2|) and l|15|l the expression for 
the variation of Casimir free energy involves the permit- 
tivity of the superconducting film e s for complex fre- 
quencies i(. Now, e s (i() cannot be given in closed form, 
and we have estimated it as follows. First, by using well 
known dispersion relations, we express e s (iQ in terms of 
the imaginary part e"(w) of the dielectric permittivity 
e a (ui) = e' s (w) + i e"(w) at real frequencies, i.e. 



e«(iC) - 1 



2 

71" ./,, 



dw 



' C 2 + u> 2 



(BIO) 



Upon using now the standard relation connecting the 
imaginariy part of the electric permittivity to the real 
part of the complex conductivity a of a metal: 



4tt 



(Bll) 



However, the sum rule holds also in the normal state of 
the film, and thus we have 



duj cr' D (uj) = 



2m 



(B15) 



where o' D {w) (see Eq. |B9j)) is the expression derived 
from Eq. (|16fl . by using Eq. I|B11|I . Upon equating 
the l.h.s. of the above two Equations, we arrive at the 
following expression for k/2: 



K 

2 



(• <_X..< 

dw (*' D (w) - *' s (w)) . (B16) 

Jo 

Upon plugging this expression for k into Eq. (|B13() , and 
by using the dispersion relation for a' L ; 



i.e. 



e D (i{) -1 = 8 



dw 



C 2 + W 2 



(B17) 



we can rewrite the r.h.s. of Eq. (|B10|) as 



e s {i() = 1 + . 



du 



C 2 + w 2 



(B12) 



We now take for a' s (w) the semianalytical BCS formula 
given in Eqs. l|B3IB8f) . Upon plugging Eq. I|B3|) into Eq. 
(IB12|I we then obtain the following formula for e s (i(): 



e a (iC) =! + ■ 



du> 



( 2 +w 2 



4k 



(B13) 



where k is given in Eq. 1|B8(I . For the purpose of a nu- 
merical evaluation of Eq. (|B13|) , it is however convenient 
to rewrite it in a different form. We recall that the co- 
efficient k is defined so as to satisfy the sum rule Eq. 
(|B2|I . Now, upon substituting the expression for cr' s (w) 
Eq. (|B3)l into Eq. l|B2)l we obtain 



+ 



dw <j' s (w) 



nne 
2m 



(B14) 



it is easy to obtain the final equation 



C 2 



dw 



w 2 (a>(w) 



H) 



C 2 +L0 2 



(B18) 

that we have used in our numerical computations. It has 
the virtue of expressing the variation in the film permit- 
tivity e s (iQ — eo(iC) across the transition, as an integral 
involving the variation of conductivity <j' s (oj) — a' D {w). 
Indeed, &' s (w) — <j' d (w) is appreciably different from zero 
only for frequencies w of order of a few times k T c /Ti, and 
goes to zero like 1/w 3 for large w's, and this ensures rapid 
convergence of the integral on the r.h.s. of Eq. IJB18J1 . 
By looking at the plots in Fig. ©, we see that, for fre- 
quencies of order 2A/h (apart from a very narrow region 
near the origin, where d s {w) has a logarithmic behavior), 
d s {w) < a n (w), and therefore it follows from Eq. (|B18(I 
that e s («C) > Cn(iC)- This relation is fully confirmed by 
our numerical computations. 
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